Neutrinoless double beta decay and pseudo-Dirac neutrino mass predictions 
through inverse seesaw mechanism 
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In the inverse seesaw extension of the standard model, supersymmetric or non-supersymmetric, while the light 
left-handed neutrinos are Majorana, the heavy right-handed neutrinos are pseudo-Dirac fermions. We show how 
one of these latter category of particles can contribute quite significantly to neutrinoless double beta decay. The 
neutrino virtuality momentum is found to play a crucial role in the non-standard contributions leading to the 
prediction of the pseudo-Dirac fermion mass in the range of 120 MeV — 500 MeV. When the Dirac neutrino 
mass matrix in the inverse seesaw formula is similar to the up-quark mass matrix, characteristic of high scale 
quark-lepton symmetric origin, the predicted branching ratios for lepton flavor violating decays are also found 
to be closer to the accessible range of ongoing experiments. 



I. INTRODUCTION: The standard gauge theory of strong, 
weak, and electromagnetic interactions has confronted numer- 
ous experimental tests while the last piece of evidence on the 
Higgs boson is currently under rigorous scrutiny at the Large 
Hadron Collider (LHC). In spite of these, neutrino oscillation 
data uncovering tiny masses of left-handed (LH) neutrinos call 
for physics beyond the standard model (SM) which is most 
simply achieved via canonical seesaw mechanism that 
requires the addition of one heavy right-handed (RH) neutrino 
per generation provided both LH and RH neutrinos are Majo- 
rana fermions 03J. Several other forms of seesaw mechanism 
J5)-01 also require Majorana fermions. Quite interestingly, on- 
going experiments on neutrinoless double beta decay (QvBfJ) 
(8D is expected to resolve the issue between Majorana y] or 
Dirac [4] nature of the neutrino In contrast to the predicted 
small contribution to the 0v/3(3 decay rate in the SM, there has 
been quite significant, or even more dominant predictions if, 
at the TeV scale, there is left-right (LR) gauge theory iflolfTlll . 
Even, attempts have been made to predict nonstandard contri- 
butions to 0^/3/3 decay rate due to the mediation of pseudo- 
Dirac neutrinos where each of them is considered to be a pair 
of Majorana neutrinos | LI , 12h . While the possibility of left- 
handed neutrinos being pseudo-Dirac has been shown to be 
highly challenging 11311 . contribution of a fourth generation 
heavy pseudo-Dirac neutrino to Qv(3f3 has been explored with 
the condition that its mass should be greater than Mz/2 lfl4ll . 
If the Dirac neutrino mass matrix occurring in seesaw formu- 
las has its left-right symmetric or quark-lepton symmetric ori- 
gin, descending from Pati-Salam symmetry [15) or SO(10) 
grand unified theory IU6I1 at high scales, then the canonical 
seesaw scale is too large to be experimentally tested by high 
energy accelerators including LHC. Alternatively, the inverse 



seesaw mechanism which requires one RH neutrino 

as well as an additional sterile fermion per generation, oper- 
ates at TeV scale and is, therefore, experimentally verifiable. 
In this framework while the LH light neutrinos are Majorana 
fermions, the RH neutrinos are pseudo-Dirac by nature having 
heavier masses. 

In this letter we show that the inverse seesaw formula ex- 
plaining the light neutrino masses and mixings permits the 
lightest of the three pseudo-Dirac neutrinos in the mass range 
(120 — 500) MeV leading to new contributions to 0i//3j3 decay 
comparable to, or much more than, those due to the exchanges 
of the light left-handed neutrinos. The neutrino virtuality mo- 
mentum dlH, \p\ ~ 190 MeV, is noted to play a crucial 
role in such new contributions. The origin of Dirac neutrino 
mass matrix is also found to be important in our estimations in 
predicting lepton flavor violating decays accessible to ongoing 
experimental searches. As our results are also applicable in 
the inverse seesaw extension of the minimal supersymmetric 
standard model (MSSM), they are consistent with gauge cou- 
pling unification at the MSSM-GUT scale, Mjj ~ 2 x 10 16 
GeV 

II. THE INVERSE SEESAW EXTENSION: As is customary to 
the implementation of inverse seesaw mechanism, we add two 
fermion singlets to each generation of the SM, with or without 
supersymmetry. While we call the first type of singlet a RH 
neutrino (Nr), the second type of singlet is named as a sterile 
neutrino (Sl) and, in the (vl, N r , Sl) basis, the 9 x 9 neutrino 
mass matrix is 111811 
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Besides the two distinct possibilities, Dirac or Majorana, very recently 
a new hypothesis has been advanced in which neutrinos could be 
schizophrenic @|. 



where Mb is the Dirac mass term of the neutrino, and M is 
the heavy Dirac mass matrix relating Nr and Sl- The matri- 
ces Md and M are in general 3x3 complex in flavor space 
whereas the /i$ is 3 x 3 complex symmetric matrix. 

Transformation from flavor to mass basis and diagonaliza- 
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tion are achieved through 

W)f = V*\v)m , 

V^M V V* = M v = Dmg{m Vi ; M c .} : 



where 



(2) 



(3) 



where \v) m = (£>i,Q) T represents the three light and six 
heavy mass states, and i and j run over the light and heavy 
mass eigenstates, respectively. With ng,Mj) <§C M, the ma- 
trix M. v can be block diagonalized to light and heavy sectors 



M D \ (M D \ 

~m) ms [la) 



M T 
M us 



(4) 



where m v has the well known inverse seesaw formula 111 811 
and Mh is the mass matrix for heavy pseudo-Dirac pairs of 
comparable masses with splitting of the order of /is- The /is 
term in the Lagrangian breaks the leptonic global symmetry, 
U(1)l, which is otherwise preserved in the standard model in 
the limit fig —> rendering all the LH neutrinos to be mass- 
less. Hence the small /is should be a natural parameter in 
this theory in the 't Hooft sense ll2lll . The above block diago- 
nalized matrices are further diagonalized through the PMNS 
matrix, U v , and a 6 x 6 unitary matrix Uh, respectively, so 
that 
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Hence, in the leading order approximation, V can be written 

as 



V 








(7) 



where X — (MdM 1 ), and all the elements in the first block 
are 3x3 matrices. 

(II. A) /is from neutrino oscillation data: The inverse see- 
saw formula in eqn. (0]l predicts light neutrino mass ma- 
trix in terms of three other matrices, Mb, M, and /is- At 
first we take Mb ~ Mi, the charged lepton mass matrix, 
which may arise if the SM originates from high scale left- 
right gauge symmetry, SU(2)l X SU(2)n x U{1)b-l X 
SU(3) C ^4 SM, where Mr » M w . Assuming the ma- 
trix M to be diagonal for the sake of simplicity and using 
M D =diag{m e ,m Al ,rn T } = {0.0005, 0.1, 1.7} GeV, we ob- 
tain /is from global fits to the neutrino oscillation data 02211 
given in TABLE [Q 



MS (GeV) 



X- 1 Mm v M T X T ] 
/ 



6.71 x 10~ 7 + 1.96 x 10 7 i -1.17 x 10~ 8 - 3.22 x 10~ 8 i -3.71 x 10" 8 - 2.03 
-1.17 x 10~ 8 - 3.22 x 10" 8 i 1.53 x 10~ 08 - 2.22 x 10~ 10 5 7 n - 1 n ~ 9 - ^ - 



x 10" 



\ -3.71 x 10~ 8 - 2.03 



x 10 



-8 



7.0 x 10~ 9 

r 



2.83 x 10 



-9 : 



7.0 X 10" 9 

-5.50 x 10~ 9 



2.83 x 10 9 i 
(-5.26 x 10" 11 




where TV = (1 — rj) U v and r\ = \XX^ is a mea- 
sure of unitarity violation. This particular structure of 
/is has been derived using, as an example, the nor- 



mal hierarchical (NH) light neutrino masses 



diag(0.00127 eV, 0.00885 eV, 0.0495 eV) and non- 
degenerate eigenvalues of M = diag {0.2, 2.6, 23.7} GeV. 
Similar analysis predicts somewhat different structures of /is 
for inverted hierarchical (IH) and quasi-degenerate (QD) pat- 
tern of the light neutrinos and can further be easily obtained 
for degenerate M\ = M% = M3 or, partially-degenerate 
Mi = Mi <C M3 after taking care of the phenomenolog- 
ical bounds |?/ ee | < 2.0 x 10~ 3 , \r}^\ < 8.0 x 10~ 4 , and 
\rj TT \ < 2.7xl0~ 3 . Our ansatz with M = diag (Mi, M 2 ,M 3 ) 



Neutrino oscillation parameters 


Globally fitted values 


Am s 2 ol [eV 2 ] 


7.58 x 1CT 5 


|Am 2 tm |[eV 2 ] 


2.35 x 10" 3 


sin 2 612 


0.320 


sin 2 #23 


0.427 


sin 2 6>i3 


0.0246 


Scp 


0.8 7v 



TABLE I: Mass squared differences, mixing angles, and CP-phase 
from global fits to neutrino oscillation data M22I1 . 
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gives 



r] = diag 



1.25 x 1CT 7 0.005 1.35 



Ml 



Ml ' Mf 



where all masses on the right hand side are in GeV. 
III. NEUTRINOLESS DOUBLE BETA DECAY PREDICTIONS 

Two separate contributions due to light and heavy neutrino 
exchanges to Oz//3/3 transition become transparent by writing 
the flavor eigenstates as linear combination of light and heavy 
mass eigenstates 

where li ~ (0, X)Uh is a 3 x 6 matrix. Then the weak charge- 
current Lagrangian can be expressed as 



£ cc = -^=W£i a ^P L v a + h.c. 
9 



^w£e a ^p L (Mam + u a j(j) + h - c - . (io) 

resulting in two different categories of Feynman amplitudes: 

• A V LL which arises from the Feynman diagram of 
Fig. Q2a) due to only light neutrino exchanges 



pseudo-Dirac neutrinos have been discussed in Sec-II. The 
half-life of 0^/3/3 transition is then found to be 



1/2 



— I^av 



V,LL 



M, 



(13) 



where /Coi> contains phase space factors plus nuclear matrix 
elements and m^ e LL (M| e LL ) represents the effective neutrino 
mass derived from light neutrino (heavy pseudo-Dirac neu- 
trino) exchanges in the mass basis. The analytic forms of the 
two effective masses have been estimated for this model as 
shown in TABLE. HH 



Effective mass 


Analytical expression 




Ml, m u% 

*>3 



TABLE II: Effective mass parameter for standard (non-standard) 
contributions due to light (heavy pseudo-Dirac) neutrino exchanges 
for 0^/3/3 decay. 



A v LL = G\Ml^, (11) 

where (p) ~ 190 MeV represents neutrino virtuality 
momentum 1 19, 20]. 

A^ LL which arises from the Feynman diagram of 
Fig-QIb) due to heavy pseudo-Dirac neutrinos, 

M c - 



A{ L = G% (W) 



ei p 2 - Mf 

^ 3 



(12) 
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-» e £ 



W 
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FIG. 1: Feynman diagrams contributing to neutrinoless double beta 
decay due to light neutrino exchanges (left-panel) and heavy pseudo- 
Dirac neutrino exchanges (right-panel). 

The mass eigenstates of heavy pseudo-Dirac neutrinos 
are , , C3"; (1 > C2 i CT) with almost degenerate pairs 
((k,Ck> k= 1,2,3) but having small mass difference /is 
between the members of the pair and the flavor states 
are (N±, N2, A3; S%, S2, S3). The mixing matrix for these 



We discuss below three different cases: 
(III. A) The standard contribution. It is well known that the 
standard contributions due to light neutrino exchanges are de- 
pendent on their allowed mass patterns; normal hierarchical 
(NH), inverted hierarchical (IH), or quasi-degenerate (QD), 



LL 



e 2la m 2 



+U, 



e3 ' 



m 3 



,LLl 



0.004 eV NH, 
0.048 eV IH, 
0.1 eV QD. 



(14) 



In our case, M e i — U e i and light neutrino exchanges in the 
mass basis gives almost the same contributions which are pre- 
sented by solid lines shown in Fig. [21 Fig- EJ and Fig. [4] 
(III. B) Mq ^> \p\: In the inverse seesaw extension under 
study, in addition to the standard effective mass parameter, the 
additional effective mass parameter for \Mq \ ^> \p\ satisfies 



M, 



<= 



This results in new contribution to 



0^/3/3 transition half-life 



T 



1/2 



= K, 
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2 / 
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\(p)\ 2 {u 



±\2 MSfcfc 



(M ( _ 
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(15) 



where fiSkk an( l ^kk are the eigenvalues of /j,s and M, re- 
spectively. One example of this case has been shown in Fig. 
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FIG. 2: Predictions of effective mass \M ee | in Of/3/3 decay with Dirac neutrino mass Md = charged lepton mass Mi and the diagonal structure 
of M for NH (IH) pattern of light LH neutrino masses as shown in the left ( right) panel. The standard contribution is shown by solid line and 
nonstandard contributions with pseudo-Dirac neutrino exchanges of different masses expressed in MeV are shown by other lines. 



|2]for — 0.5 GeV where the predicted effective mass pa- 
rameter is nearly 3/2 (4) times larger than the standard predic- 
tion for NH (IH) case. 

(III. C) ■ ~ \p\: In this region where different allowed val- 
ues of Mf are of the order of neutrino virtuality momentum 
\p\ ~ 190 MeV, the new contribution to neutrinoless double 
beta decay due to heavy pseudo-Dirac neutrino exchange is 
found to be more dominant than the standard contribution and 
the Oz^/3/3 transition half-life is given below 



1/2 



pscudo — Dira(^ 




The predicted new values of the effective mass parameters 
arising solely due to pseudo-Dirac neutrino exchanges have 
been shown in Fig. [2]in the left-pannel (right-pannel) for NH 
(IH) patterns of the light neutrino masses, respectively, where 
M$ x = (0.15 - 0.5) GeV. It is quite clear from the plots that 
even for — 0.25 GeV or, 0.5 GeV, the new contribu- 
tions are 3-6 times larger than the standard ones. While for 
the value of — 0.18 GeV, the contribution is nearly 100 
times larger shown in Fig. [2] This large enhancement occurs 
as approaches the vicinity of the neutrino virtuality mo- 
mentum, \p\ ~ 190 MeV. We point out that such important 
effects of pseudo-Dirac neutrino masses are found for the first 
time in this work. 

IV. LEPTON FLAVOR VIOLATION WITH DOMINANT Oi//3/3 
DECAY RATE: We have clearly shown that the predicted non- 
standard contributions to neutrinoless double beta decay rate 
are dominant for the lightest allowed pseudo-Dirac neutrino 
mass M < ^ 1 ~ (0.15 — 0.5) GeV. However, because of the diag- 
onal nature of Md and assumed structure of M, the branching 
ratios for lepton flavor violating (LFV) decays, /j, —> e + 7, 
t — > e + 7, and r — > fi + 7 are as small as the SM predictions. 



In the next two examples we adopt plausible parametrization 
predicting significantly larger contribution to these branch- 
ing ratios while retaining the dominant contributions to 0i/j3(3 
transition. 

(IV. A) M D ~ Mi with non-diagonal M : We generate non- 
diagonal matrix M to satisfy the existing phenomenological 
bounds on the elements of 77 [23], 



|77 ee | < 2.0 x 10~ 3 , < 3.5 x 10" 

|Tfer| <8.0X 10~ 3 , < 8.0 X 10" 



< 5.1 x 10" 



['h 



< 2.7 x 10" 



(17) 



Using the parametrization of the type used in ref. [24], M can 
be expressed as 



M 



1 



O t M d 



1 1 r 



M 



1 



O t M d 



= U =V^V 



(18) 



where O is the matrix diagonalizing \rj\ and V is an arbitrary 
unitary matrix. Choosing, for the sake of simplicity, V = I3 
and we note that that a lightest pair with Mf . ~ 0.16 GeV, in 
the vicinity of neutrino virtuality momentum, can be achieved 
by suitable rescaling, e.g. ri a p —> r/ Q| g/(1500). After this 
scaling we find 



M(GeV) = 



1 



2-q 

/0.092i 
0.217 
^0.074 



O t M 



D 

14.08? 
-70.36 
8.866 



(19) 



-383.7^ 
-80.39 
320.5 



where 



O 



1 0.5874i 0.5446 0.5987 
0.4284i -0.8368 0.3409 
y-0.6866i -0.0562 0.7248 
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FIG. 3: Same as Fig.[2]but now with non-diagonal structure of M and reduced values of nonunitarity matrices rj as described in the text. 



With the allowed mass eigenvalues for the heavy pseudo- 
Dirac neutrinos, M ( = diag{0.159, 72.0, 506.4} GeV, the 
predicted branching ratios for lepton flavor violating decays 
areJH 



Br(/w e + 7 ) = 1.56 x 10~ 26 , 
Br (r ->• e + 7) = 5.79 x 10" 27 , 
Br(r^ M + 7) = 1-10 x 10~ 18 . 



Although all the three branching ratios are much smaller than 
their corresponding experimental upper limits H27L they are 
considerably larger than the SM predictions. However we 
note below that with Mo similar to M u , the up-quark mass 
matrix, a phenomenon underlying the possible origin of SM 
from Pati-Salam lfl5ll or SO(10) model, LFV decays have 
much larger predicted values, accessible to ongoing experi- 
mental searches, while similar predictions on dominant Of/5/3 
decay are maintained. 

(IV. B) M D ~ M u and GUT connection: In this case 
Dirac neutrino mass matrix is approximated to be up-quark 
mass matrix, which originates if the high scale symme- 
try is Pati-Salam or SO(10) GUT, SU(2) L x SU(2) R x 



SU(4) C or SO(10) ^ SM. Using running masses 
(m u ,m c ,m t ) = (0.00233, 1.275, 160) GeV and Cabbibo- 
Kobayashi-Maskawa mixing matrix, Vckm [28], 



M D (GeV) 

/0.067-0.004i 
0.302- 0.022z 
y 0.55 - 0.53?; 



M u — VckmM u vJ km 



0.302- 0.022?; 0.55 - 0.53?; 

1.48 -O.Oi 6.534 - 0.00H I .(21) 
6.534- 0.0009?; 159.72 + O.Oi 



At first using the phenomenological bounds from eqn. (fTTt 
and saturating our ansatz for rf a p = \r)\ max /n, n = 1, ..5 , we 
search for matrix M through eqn. (fT8l which gives /is from 



eqn. ©. We obtain for n = 4 

M(GeV) 

-5.9-3.45i 0.2 - 60.72?; 

-10.44 + 2.13* -60.9 - 0.5i 

7.08- 5.09?; 70.86 - 0.18i 



5.15 - 1760?; 
-598.0- 12. 16t 
1547- 4.15i 



(20) ^ s (eV) 



/-3.42-0.5H 
-1.92 + 5. 55i 
v 0.28 - 1.92i 



-1.92 
39.1 4 
-12.1 



-5.55i 0.28- 
5.68z -12.1 
-O.lli 4.10- 



1.92i 
-O.lli 
0.68z 



(22) 



•(23) 



Our predictions on numerical values of the effective mass 
parameter for 01/(3/3 are shown in Fig. |4]for NH, IH and QD 
cases. For NH light neutrinos we find that the predicted value 
of \M ee \ is increased by a factor 3 for 77 = \r)\ max /3, corre- 
sponding to lightest pair M^ l = 131 MeV, while the incre- 
ment is 10 times for rj = \r\\ max j\ and Mf, = 152 MeV, 
and 30 times for 77 = \r]\ max /5 and M Cl = 169 MeV We 
find that the enhancement survives as long as lightest pair 
M (l ~ 120 - 350 MeV. For the IH light neutrino masses the 
results are similar as shown on the right panel of Fig. |4] The 
branching ratios for lepton flavor violating decays predicted 
in this scenario with M c = (0.152, 39.5, 2426) GeV are 



e + 7 ) = 3.6 x 10~ 13 : 
e + 7 ) =4.2 x 10~ 14 : 



Br(/x 
Br(r 

Br(r->-M + 7) = 3.3 x 10~ 12 



(24) 



while the present experimental limits at 90% C.L. on these 
branching ratios are Br (fi — > e + 7) < 1.2 x 10~ n , 
Br (r -> e + 7) < 3.3 x 10" 8 , and Br (r -> \i + 7) < 4.4 x 
10 -8 B27I1 . The projected reach of sensitivity in the future 
is Br (r — > e + 7) , Br (r — > j jl + 7) < 10 -9 and specifically 
Br (/i -> e + 7) < 10" 19 0. 

The predicted nonstandard contributions to 0v(3/3 transition 
are shown in the left-panel for NH and in the it right-panel 
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FIG. 4: The effective mass parameter, |M ee |, prediction for Ou/3/3 decay due to light (Solid line) and pseudo-Dirac (dashed lines) neutrino 
exchange where Mr> — M u and rj — \r) max \/n, n — 1, 3, 4, 5 as discussed in the text. 



for IH case of Fig. [4] In view of the -dependent en- 
hancements of 0^/3/3 decay rates discussed above it is tempt- 
ing to search for the possibility of the lightest pseudo-Dirac 
neutrino mass which we perform by the replacement M£ — > 
M| + iM^Tx, where T\ corresponds to plausible value of 
width of the particle. Using, for example, Ti ~ 0.1 keV, our 
predictions are presented by solid curve in Fig. |5]for NH light 
neutrino masses where the resonant behavior is clearly exhib- 
ited around M Cl = 190 MeV. 




170. 180. 190. 

Mr 



200. 210. 
(MeV) 



220. 230. 



FIG. 5: Variation of effective mass |M ee | prediction as a function 
of lightest pseudo-Dirac neutrino mass M^ t where us matrix has 
been determined for NH light neutrino masses. The resonance peak 
is at M( 1 ~ 190 MeV as shown by the solid line. For comparison, 
the standard contribution with NH masses is shown by the dashed 
horizontal line. 

V. DISCUSSIONS AND CONCLUSION: In this letter we have 
shown that in the inverse seesaw framework of the standard 
model, the lightest of the pseudo-Dirac neutrino could be of 
0(100) MeV in concordance with tiny left handed neutrino 
masses and the oscillation data. This pseudo-Dirac neutrino 
mass being in the vicinity of the neutrino virtuality momen- 
tum \p\ ~ 190 MeV, gives very significant non-standard con- 



tributions to 0v(3f3 decay rates, even far exceeding the stan- 
dard contributions. The Dirac neutrino mass possibly origi- 
nating from high scale Pati-Salam symmetry or SO(10) grand 
unification, plays a crucial role in determining dominant con- 
tributions to Qi//3/3 decay rates simultaneously with LFV de- 
cays with predicted branching ratios accessible to on going 
search experiments. The underlying mechanism provides 
three distinct platforms for its falsifiability (i) 0vf3f3 decay 
rates, (ii) determination of light pseudo-Dirac neutrino mass 
M C:1 ~ 120 — 500 MeV, and (iii) the three predicted branching 
ratios of eqn. (124) . As all our results are applicable in the case 
of inverse seesaw extended supersymmetric standard model, 
they are also consistent with gauge coupling unification at the 
MSSM-GUT scale, M v ~ 2 x 10 16 GeV. The Pati-Salam or 
SO(10) completion of the model discussed in Sec. IV.B will 
be reported elsewhere in future publication ()26|. 
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